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Abstract. Stress enhancement in the vicinity of brittle cracks makes the 
macro-scale failure properties extremely sensitive to the micro-scale material 
disorder. Therefore: (i) Fracturing systems often display a jerky dynamics, so- 
called crackling noise, with seemingly random sudden energy release spanning 
over a broad range of scales, reminiscent of earthquakes; (ii) Fracture surfaces 
exhibit roughness at scales much larger than that of material micro-structure. 
Here, I provide a critical review of experiments and simulations performed in this 
context, highlighting the existence of universal scaling features, independent of 
both the material and the loading conditions, reminiscent of critical phenomena. 
I finally discuss recent stochastic descriptions of crack growth in brittle disordered 
media that seem to capture qualitatively - and sometimes quantitatively - these 
scaling features. 
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1. Introduction 

Understanding how materials break is of fundamental and technological interest. As 
such, it has attracted much attention from engineers and physicists. Since the pioneer 
work of Griffith [1], a coherent theoretical framework, the Linear Elastic Fracture 
Mechanics (LEFM) was developed. It states that, in an elastic material, crack initiation 
occurs when the mechanical energy released by crack advance is sufficient to balance 
that needed to create new surfaces [1]. Then, once the crack starts to grow, inertial 
effects should be included in the energy balance (see e.g. [2]). This continuum theory 
has proven to be extremely powerful to describe ideal homogeneous media, but fails 
to capture some features observed in heterogeneous ones. 

While the mechanical energy released during crack growth is well determined 
by continuum theory, the dissipation processes occur within a tiny zone at the 
crack tip. This separation in length-scales, intrinsic to fracture mechanics, makes the 
failure properties observed at the continuum length-scale extremely sensitive to the 
microstructural disorder of the material. Consequences include important fluctuations 
in the strength exhibited by different samples of the same material [3] , size effects in 
the strength [4, 5], crackling dynamics with seemingly random violent events of energy 
release prior and during the failure [6] and rough crack paths [7], among others. 

I will focus in this review on two of these aspects, namely the crackling dynamics 
evidenced in brittle failure and the morphology of fracture surfaces. Experimental 
observations performed in this context report the existence of scale invariant laws 
[8, 6, 9, 10]. This has attracted much attention from the statistical physics community 
over the last 25 years [11, 12, 13]. In this respect, lattice models such as Random 
Fuse Models (RFM) [14] were shown to capture qualitatively these scaling features in 
minimal models, keeping only the two key ingredients responsible for the complexity 
in material breakdown: The material local disorder and the long-range elastic load 
redistribution after a local failure event. More recently, it was proposed to extend 
LEFM descriptions to heterogeneous brittle materials by taking into account the 
microstructural disorder through a stochastic term [15, 16]. As a result, the onset 
of crack propagation is analogue to a dynamic phase transition between a stable 
phase, where the crack remains pinned by material heterogeneities, and a moving 
phase, where the mechanical energy available at crack tip is sufficient to make the 
front propagate. As such, it exhibits universal - and predictable - scaling laws that 
reproduce fairly well the intermittency statistics observed in crack dynamics [ ] and 
the morphological scaling features exhibited by fracture surfaces [18]. 

This review is organized as follow. First, I provide a brief introduction to standard 
LEFM theory in section 2 and discuss some of its predictions in term of fracture 
dynamics and crack roughness. Then, I review in section 3 the experiments and fields 
observations that have allowed to characterize the intermittent dynamics observed in 
many fracturing systems, from earthquakes to laboratory fracture experiments. Section 
4 is devoted to a critical overview of the morphological scaling features evidenced 
experimentally in various materials. RFM models and their predictions in term of 
crackling dynamics and crack roughness are discussed in section 5. Section 6 discusses 
the stochastic descriptions of crack growth and highlights their predictions in term 
of scale-free distributions, scaling laws and scale-invariant morphology for fracture in 
perfectly brittle disordered materials. Finally, section 7 outlines the current challenges 
and possible perspectives. 
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2. Continuum theory of crack growth 

2.1. Stress concentration 

Crack propagation is the basic mechanism leading to the failure of brittle materials. 
The crack motion in a brittle homogeneous material is classically analyzed using 
methods from elastodynamics. Let us first consider the case of a straight running 
crack embedded in a linear elastic material under tension. As first noted by Orowan 
[19] and Irwin [20], the stress field is singular in the vicinity of the crack tip and 
can be written as (see e.g. [2]): 



where v is the crack speed, (r, 9) are the polar coordinates in the frame (e x , e y ) centered 
at the crack tip (see Figure la), and Fij(9,v) and Gij(9,v) axe some non-dimensional 
universal functions indicating the angular variation of the stress field, and its variation 
with v. The two prefactors Kj and T, called Stress Intensity Factor and T -stress 
respectively, only depend on the applied loading and the specimen geometry. They 
determine the intensity of the singular and non-singular parts of the local field. 

This can be generalized to any loading conditions. It is convenient to decompose 
them into the sum of three modes (Figure lb): 

• Mode / (tensile mode) corresponds to normal separation of the crack under tensile 
stresses. 

• Mode II (shearing mode) corresponds to a shear parallel to the direction of crack 
propagation. 

• Mode III (tearing mode) corresponds to a shear parallel to the crack front. 

In the vicinity of the crack tip, the stress field is then written as a sum of three terms 
the shape of which is given by the right hand term of equation 1 with prefactors 
{Kj,Tj}, {Kjj,Tjj} and {Kjjj,Tjjj} associated to the tensile, shearing and tearing 
modes, respectively. As we shall see in section 2.3, a crack propagating in an isotropic 
solid chooses its orientation so that it makes shear vanish at its tip. As a consequence, 
in most of the crack problems discussed thereafter, we will consider solids loaded 
dominantly in tension where the mode II and III perturbations will be subordinated 
to that in mode I. Two exceptions will be discussed later in this review, namely 
earthquakes problems and rocks broken under compression. In both cases shear and/or 
tear fracture modes are dominant. 

2.2. Crack propagation: Griffith theory 

Because of the stress singularity, the mechanical energy released as fracture occurs is 
entirely dissipated within a small zone at the crack tip, referred to as the Fracture 
Process Zone (FPZ). In nominally brittle materials, the absence of an outer plastic 
zone can be assumed. According to Griffith's theory [1], the onset of fracture is then 
reached when the amount of elastic energy released by the solid as the crack propagates 
by a unit length is equal to the fracture energy T defined as the energy dissipated in 
the FPZ to create two new fracture surfaces of unit length. The form taken by the 
stress field at the crack tip (Equation 1) allows relating the mechanical energy release 
G at the onset of crack propagation (v — 0) to Kf. 



F ij (6,v)+TG ij (d,v) + 0(Vt), 



(1) 




G = 



Kj 



(2) 



E ' 
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Figure 1. (a) Sketch and notations describing the stress field in the vicinity of a 
slit crack tip in a two-dimensional medium under tension, (b) The three modes 
of fracture: Tensile mode (Mode I), shearing mode (Mode II) and tearing mode 
(mode III). 



where E is the Young modulus of the material. The Griffith criterion for fracture onset 
then reads: 

G > T or equivalently Ki > K c with K c = \ r ET, (3) 

where the so-called toughness K c is, like T, a material property to be determined 
experimentally. Then, once the crack starts to propagate, its motion is governed by 
the balance between the mechanical energy flux into the FPZ and the dissipation rate 
Tv. From the form taken by the stress field at the crack tip (Equation f), the energy 
flux can be related to Ki and v, which yields (see e.g. [2]): 

A(v)fy=T with A(v)~(l-— V (4) 

E V c rJ 

where cr refers to the Rayleigh wave speed in the material. The dynamic fracture 
regime reached when v is on the order of Cr is beyond the scope of this review. An 
interested reader can refer to the reviews of Ravi-Chandar (f998) [2f] and Fineberg 
and Marder (1999) [22] for a detailed presentation of this regime. In the sequel, I will 
focus on the slow crack growth regime where v <C cr. Then, the equation of motion 
can be rewritten as: 

—v~Kj — K c or equivalently — v ~ G — T, (5) 

where the effective mobilities [i and [t! are given by fi — 2cr/K c and // = cr/T, 
respectively. 

2.3. Crack path: Principle of Local Symmetry 

Finally, to complete the continuum theory of crack growth in ideal homogeneous brittle 
materials, one has to add a path criterion. This is provided by the Principle of Local 
Symmetry (PLS) of Goldstein and Salganik (1974) \2?>\ which states that a moving 
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crack progresses along a direction so as to remain in pure tension. In 2D systems, the 
crack is loaded by a combination of mode I and II only and the PLS implies that the 
direction of crack propagation is chosen so that Ku = 0. In 3D systems, the crack 
load can also contain a mode III component. Then, to cancel Km and to propagate 
in pure mode I, the crack front would have to twist abruptly around the direction of 
propagation, which would yield unphysical discontinuities in the crack path. In this 
situation, the front is commonly observed to split into many pieces and to form 'lances' 
[24] the origin of which remains presently highly debated (see e.g. [25, 26] for recent 
works on this topic). 

2-4- Continuum mechanics predictions 

There are numerous consequences of the Linear Elastic Fracture Mechanics (LEFM) 
summarized above. Here I will focus on two of them: (i) The equation of motion 
(Equation 5) predicts regular and continuous dynamics for crack propagation; (ii) 
The PLS predicts smooth surfaces at continuum scales, i.e. at scales over which the 
mechanical properties of the considered brittle material are homogeneous. As one will 
see in the two next sections, there are numerous experimental observations that have 
been accumulated over the past three decades that contradict both of these predictions. 

3. Crackling dynamics in fracture 

3.1. Earthquake statistics 

Contrary to that predicted by LEFM, crack growth in heterogeneous brittle materials 
often displays a jerky dynamics, with seemingly random discrete jumps of a variety of 
sizes. The most obvious signature of this crackling dynamics can be found in the seismic 
activity of faults: Earth responds to the slow shear strains imposed by the continental 
drift through a series of violent impulses, earthquakes, spanning over many orders of 
magnitude, from barely noticeable to catastrophic ones (Figure 2a). The distribution 
of radiated seismic energy presents the particularity to form a power-law (Figure 2b) 
with no characteristic scale: 

P{E) oc E~ p , (6) 

where the exponent (3 ~ 1.7 was proven to be very robust over different regions [27, 28]. 

Earthquake sizes are most commonly quantified by their magnitude M. The first 
(and most popular) magnitude scale is the one introduced by Richter (1935) [29] 
that relates M to the logarithm of the maximum amplitude measured on a Wood- 
Anderson torsion seismometer at a given distance from the earthquake epicenter. In 
this respect, M is relatively easy to measure and it provides information that is directly 
useful for various engineering applications. Magnitude and radiated seismic energy of 
earthquakes are usually related via the Gutenberg- Richter empirical relation [30]: 

log 10 J5 = 1.5M + 11.8, (7) 

where E is expressed in Joule. While this relation is considered to be reasonably 
accurate for earthquakes of small and intermediate sizes, its validity is questioned for 
great earthquakes [31]. Using this energy-magnitude relation in equation 6 leads to 
the Gutenberg-Richter frequency-magnitude statistics: 



p(M) oc exp(-bM), 



(8) 
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where b and j3 are related through: /3 = 6/1.5 + 1. 

The b measurement has been one of the most frequently discussed topics in 
seismicity studies (see e.g. [27] for a recent review). Observed b values of regional 
seismicity typically fall in the range 0.7 — 1.3 and usually take a value around 1 [27]. 
Its universality (or not) remains an open question [27], made particularly tough since 
b measurement can be affected by various factors as e.g. the type of magnitude scale 
chosen, the completeness and the magnitude range of the collected data, to name a 
few. 

The radiated seismic energy is also found to scale with the rupture area S of 
earthquakes as [32, 33]: 

EocS 1 - 5 . (9) 

Using this relation in equation 6 yields the following distribution for rupture area: 

p(S) oc S~ T , (10) 

with t = 1.5/3 - 0.5 ~ 2. 

The time occurrence of earthquakes exhibits also scale-free features expressed 
by Omori's law [8], which states that a main earthquake is followed by a sequence 
of aftershocks the frequency of which decays with time as t~ a , with a ~ 1. More 
recently, driven by the availability of complete seismic catalogs (see e.g. [34]) gathering 
the occurrence time, the magnitude, and the 3D location of earthquakes, many 
studies [35, 36, 37, 38, 39, 40, 41] were performed to characterize the spatio-temporal 
organization of earthquake events. In particular, Bak et al (2002) [ ] showed that 
the distribution of waiting times P{At) between earthquakes of energy larger than E 
within a zone of size L obeys to the following generalized expression of Omori's law: 

P E ,L( At ) = At -Q 7(At/Ai ) with At ocl/(E- p L d f), (11) 

where f3 is the exponent defined in Eq. 6, a ~ 1 is the Omori exponent, df ~ 1.2 is 
a fractal dimension characterizing the spatial distribution of epicenters and f(x) is a 
function constant for x <C 1, and decaying rapidly with x for More recently, 

Davidsen and Paczuski (2005) showed [39] that the spatial distance Ar between the 
epicenters of two successive earthquakes of energy larger than E within a zone of size 
L exhibits scale-free statistics: 

P E , L (Ar) = (Ar/L)- 5 /(Ar/L), (12) 

where 6 ~ 0.6 and f{x) is a function constant for x -C 1, and decaying rapidly with x 
for i>1. 

Since these scale-free distributions are observed universally on Earth, 
independently of the considered area of investigation and the period over which the 
analysis was performed, it was suggested [42] that similar statistical features should 
be reproducible in laboratory experiments. In this context, time series of Acoustic 
Emission (AE) were recorded in various rocks loaded under uniaxial compression up 
to (shear) fracture [42, 43, 44, 45]. As a result, the waiting time distribution P(At) 
was found [ ] to obey Omori's law given by equation 11 with an exponent a ~ 1 very 
close to that observed in earthquakes. 
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Figure 2. (a) Energy radiated by earthquakes of magnitude M > 3 for each day 
of the year 2005, in a region of California spanning 32°N — 37°N latitude and 
114°W — 122°W longitude, (b) Distribution in energy of for earthquakes having 
occurred in this region from 1987 to present. The axes are logarithmic. The dash 
line is the power-law given by equation 6 with fi = 1.7. (c) Distribution of waiting 
time At between two successive earthquakes of magnitude larger than 2.5 within a 
zone of size L = 1° in this region from 1987 to present. The axes are logarithmic. 
The dash line is the Omori's law (Equation 11) with a = 0.9. (c) Distribution of 
jumps size Ar between two successive earthquakes of magnitude larger than 2.5 
within a zone of size L = 1° in this region from 1987 to present. The axes are 
logarithmic. The dash line is the law of Davidsen and Paczuski (Equation 12) with 
8 = 0.7. Data were taken from seismic catalogs available at Southern California 
Earthquake Data Center (SCEDC) [34]. 



3.2. Acoustic emissions in tensile failure 

Signature of crackling dynamics is also evidenced in the AE that goes with the 
mode I failure of many brittle materials: The distributions of energy and silent time 
between two successive events have been computed in many fracture experiments 
[46, 47, 48, 45, 49, 50, 51] and found to follow power-laws similar to what observed 
in earthquakes (See equations 6 and 11). It should be emphasized however that the 
relation between AE energy and released elastic energy remains largely unknown (see 
e.g. [52] for recent work in this context). The exponents (3ae and a associated with 
AE energy and silent time were reported to depend on the considered material: They 
were found to be e.g. {/3ae — 1-5, a ~ 1.9} in wood [46, 47], {Pae — 2,a ~ 2.7} 
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in fiberglass [46, 47], {Pae — 1-3 — 1.8, a ~ 1 — 1.5} in polymeric foams [50, 51]|, 
{(3ae — 1, a — 1.25} in sheets of paper [48]. As a synthesis, the values for 0ae and a 
measured in mode I fracture experiments range typically between 1 and 2. Santucci et 
al (2004) [53, 54] have also observed directly the slow sub-critical intermittent crack 
growth in 2D sheets of paper, up to their final breakdown. In these experiments, the 
crack tip is found to progress through discrete jump the size of which is power-law 
distributed with an exponent close to 1.5 up to a stress-dependent cut-off that diverges 
as a power-law at K c . The waiting time At between two jumps is found to exhibit 
power-law distribution with Omori's exponent a close to 0.7 for At smaller than the 
mean value, and a = 2 for At larger than the mean value. 

It is worth to notice that most of AE fracture experiments reported in the 
literature are non-stationary: Usually, one starts with an intact specimen and loads 
it up to its catastrophic failure. In these tests, the recorded AE reflects more the 
micro-fracturing processes preceding the initiation of the macroscopic crack than the 
growth of this latter. More recently, Salminen, Koivisto and co-workers (2006) [55, 49] 
investigated the AE statistics in a steady regime of crack propagation in experiments 
of paper peeling. The distribution of AE energy follows a power-law with an exponent 
Pae — 1-8 ± 0.2 significantly higher than that observed in standard tensile tests 
starting from an initially intact paper sheet [48]. Nevertheless, the silent time between 
successive events was found to display statistics and correlations significantly more 
complex than that of Omori'law. 

3.3. Intermittent crack growth along planar heterogeneous interfaces 

AE methods are very powerful to localize micro-fracturing events in time. On the 
other hand, to relate quantitatively the AE energy to the local crack dynamics or 
the mechanical energy released through these events remains a rather difficult task 
(see e.g. [52] for recent work in this context). This has motivated Maloy, Santucci 
and co-workers to investigate experimentally the dynamics of crack propagation in a 
simpler configuration, namely the one of a planar crack propagating along the weak 
heterogeneous interface between two sealed transparent Plexiglas plates [56, 57]. Using 
a fast camera, they observe directly crack propagation (Figure 3a, from [57]). 

To characterize the local dynamics, they computed the time spent by the crack 
front at each location (z, x) of the recorded region. A typical gray-scale image of this so- 
called waiting-time map, w(z,x) 1 is shown in Figure 3b (from [57]). The intermittency 
reflects in the numerous and various regions of gray levels. The local speed v(z,x) of 
the crack front as it passes through a particular location (z, x) is then proportional 
to l/w(z,x). Quakes are then defined [57] as connected zones where v{z,x) > c{v) 
where () denotes averaging over all locations and c is a constant ranging from 2 
to 20 (Figure 3c, from [57]). Recent analyses performed by Grob et al (2009) [58] 
reveal many similarities between the statistics of these experimental quakes and that 
of real earthquakes: First, the distribution of quake area, S w , follows a power-law, 
characterized by an exponent t w ~ 1.7 [57], which was suggested [58] to be analogous 
to the energy distribution characterized by an exponent /3 ~ 1.7 observed in real 
earthquakes (Equation 6); Second, the distribution of time occurrence and spatial 
distance between the epicenters of two successive events were found to be the same in 
both experimental quakes and real earthquakes (Equations 11 and 12, respectively). 

% In polymeric foams, both /3ae and a were also found to depend on the loading conditions (load- 
imposed or displacement-imposed) and on the temperature [50, 51] 
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Figure 3. (a) Typical image recorded by Mal0y et al of the crack front 
propagating along the weak heterogeneous interface between two sealed Plexiglas 
plates. The solid line represents the interface separating the uncracked (in black) 
and cracked (in gray) parts extracted after image analysis, (b) Gray scale map 
of the waiting time matrix. The darker parts show the longer waiting times, (c) 
Spatial distribution of quakes (in white) corresponding to clusters with waiting 
time 10 times smaller than the mean one. (Courtesy of K.-J. Mal0y, from [57]). 



All these distributions were found to be independent of the clipping parameter c used 
to define the quakes in the waiting-time maps [57, 58]. 

Let us draw some main conclusions from the experimental and fields' observations 
reviewed here: 

(i) Failure of brittle heterogeneous materials displays crackling dynamics, with 
sudden random events of energy release and/or jumps in the crack growth. The 
event statistics is characterized by various power-law distributions in e.g. energy, 
size or silent time between successive events. This crackling dynamics is observed 
in the micro-fracturing processes occurring prior to the initiation of a macroscopic 
crack as well as during the slow propagation of this latter. 

(ii) The exponents associated with the distribution of energy and silent times between 
the micro-fracturing events preceding the initiation of a macroscopic crack are 
generally observed to depend on the material, loading conditions, environment 
parameters... 

(iii) The intermittency observed during slow crack propagation in an heterogeneous 
brittle medium loaded under tension shows strong similarities with the earthquake 
dynamics of faults. The exponents associated with the different power-law 
distributions seem to take universal values. 
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Figure 4. Distribution P(S W ) of the area of the quakes detected on the waiting 
time maps computed in the model experiment of Mal0y et al (see Figure 3). The 
notation S = S w was the one used in ref. [57]. The various symbols correspond to 
various values of mean front velocity (v) and pixel resolution a. Inset: Distribution 
of quake area for a wide range of clipping parameter c. The axes are logarithmic. In 
both plots, the straight line corresponds to a power-law distribution with t w ~ 1.7. 
(Courtesy of K.-J. Mal0y, from [57]). 



Finally, it is worth to mention that crackling dynamics has also been observed, at 
the nanoscale, in cleavage experiments performed on pure crystals like mica [V ] and 
sapphire [60, 61]. In the second case, both the distributions of released energy and 
silent time between two successive nanofracture events have been computed. They are 
found to obey power-law similar to that observed in AE during the fracture of brittle 
heterogeneous materials or in earthquakes. 

4. morphology of cracks 

We turn now to crack roughness. The morphology of fracture surfaces has been widely 
investigated over the last century - fractography is now routinely used to determine the 
causes of failure in engineering structures (see e.g. [7] for a recent review). In section 
2.3, we saw that LEFM states that crack propagation in brittle materials remains in 
pure tension as long as v remains small. This would lead to smooth fracture surfaces. 
Experimentally, the failure of disordered media leads to rough surfaces at length scales 
much larger than that of the microstructure (Figure 5). 

4--1- Self-affine scaling features of crack profiles 

Starting from the pioneer work of Mandelbrot et al (1984) [9], many experiments 
have revealed that fracture surfaces exhibit self-affine morphological scaling features. 
In other words, crack profiles such as the one shown in Figure 6a are statistically 
invariant through the affine transformation (r,h) — > (Ar, X^h) where r and h(r) refer 
to the in-plane and out-of-plane coordinates respectively. The exponent £ is called the 
Hurst or roughness exponent. Experimentally, this affine invariance can be checked 
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Figure 5. Typical fracture surfaces in (a) quasi-crystals (STM, Courtesy of P. 
Ebert, from [62]), (b) glass (AFM, from [63]), (c) mortar (Contact profilometry, 
courtesy of S. Morel, from [64]) and (d) glassy ceramics made of sintered glass 
beads (Interferometric profilometry). In all these materials, the roughness extends 
over length-scales much larger than that of the microstructure. 



by computing the distribution P& r (Ah) of height increments between two points 
separated by a distance Ar (Figure 6b): For self-affine crack profiles, Pa?- (Aft.) takes 
the following form (Figure 6c): 

PAr(Aft) = Ar- c f(Ah/Ar c ) (13) 

The shape of the generic distribution / is often found to be very close to a 
Gaussian (see e.g. [65, 66] and Figure 6c). Then, £ allows one to characterize entirely 
the statistics of crack roughness. In this respect, its determination was an important 
issue and lead to many studies over the last 25 years. Besides the full computation of 
the distribution P/^ r (Ah) that presents the drawback to necessitate a large amount of 
data, several other methods were developed to measure (. They will not be reviewed 
here. An interested reader may refer to the work of Schmittbuhl (1995) [67] for more 
information on these methods and their limitations. In the following, we will use one 
of the most common ones, which consists in computing the height-height correlation 
function: 



Ah(Ar) = ({h{r + Ar) - h{r)) 2 )l /2 , 



(14) 
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where (} r denotes average over all positions r along a given profile. For a self-affine 
profile, one gets: 

Ah(Ar) oc Ar c (15) 

It is important to emphasize that this method, - as well as many others -, allows to 
measure the roughness exponent provided that the profile is actually self-affine, but do 
not ensure this self-affinity by itself. 

Possible universality in the morphological scaling features of cracks was first 
mentioned by Termonia and Meakin (1986) [68]: They made use of a minimal 2D 
molecular model of material failure to address the question and observed that the 
fractal dimension of the final crack was independent of the elastic constants over a 
wide range of values. Bouchaud et al (1990) [69] were the first to suggest a universal 
value for the roughness exponent £ on the basis of experimental observations: In their 
seminal series of fractography experiments, the authors observed that £ ~ 0.8 in 
four specimens of aluminium alloy on which different heat treatments had conferred 
different microstructures and toughness. This universality seemed to be corroborated 
by Maloy et al (1992) [70] who observed similar ( values in six different brittle 
materials. Thereafter, many experimental works were focused on the measurement 
of £ - the review written by Bouchaud (1997) [ ] synthesizes most of them up to 
1997. Reported £ values are between 0.5 and 0.9 in metallic alloys [9, 69, 71, 72], 
around 0.9 in graphite [70], around 0.75 in porcelain [70], around 0.85 — 0.9 in Bakelite 
[70], between 0.5 and 0.8 in oxide glasses [73, 74], around 0.8 — 0.9 in wood [ ] and 
in ice [76], around 0.45 — 0.5 in Sandstone [77, 78, 66], around 0.45 — 0.5 in Basalt 
[78], around 0.7 — 0.8 in mortar [79], and around 0.4 — 0.5 in glassy ceramics made 
of sintered glass beads with various porosities [80, 81]. As a result, these fractography 
experiments performed over a period spreading from 1984 to 2006 leaded to mitigated 
opinions on the relevance or not of this universality [71, 82]. As we will see in the two 
next sections, some recent observations allow to reconciliate this apparent spreading 
in the measurements of ( with the scenario of universal morphological scaling features. 

All the measurements reported in the preceding paragraph concern the roughness 
of fracture surfaces resulting from the failure of 3D specimens. Crack morphology was 
also explored in 2D geometries. Crack lines in paper broken in tension were also found 
to be self-affine with a roughness exponent ( ~ 0.7 [83, 84, 65]. Note that a small 
but clear difference was found in ( between the slow (subcritical) and the fast growth 
regime [85]. Crack front in the interfacial fracture experiments presented in section 
3.3 were also reported to exhibit self-affine scaling features. The associated roughness 
exponent was found to be ~ 0.6 [86, 87]. 

Let us end this section by mentioning that this picture of simple self-affine 
cracks was recently questioned. First, multi-scaling, characterized by non-constant 
scaling exponents C, q between the higher order height-height correlation function 

Ah q (Ar) = ((h(r + Ar) — h(r)) q )V q and Ar, was invoked in both 2D geometries 
[88], and 3D geometries [89]. This multi-scaling seems however to disappear at large 
scales [84, 90, 65]. Second, fracture surfaces exhibit anomalous scaling [91, 75]: The 
introduction of an additional global roughness exponent, Cgiobi is necessary to describe 
the scaling between the global crack width and the specimen size. Experiments reveal 
that (giob depends on the material and on the fracture geometry [91, 75]. 




Figure 6. (a) Typical fracture profile as observed via AFM in silica glass broken 
under stress corrosion. This profile was taken perpendicular to the direction of 
crack growth, (b) Distribution of the increments Ah between two points of such 
profiles for various values of Ar. (c) Data collapse of this distribution using Eq. 13 
with £ ~ 0.8 which shows that the profile is self-affine with a roughness exponent 
C = 0.8. 



4-.2. Family- Vis eck scaling of fracture surfaces 

To uncover the primary cause of the apparent spreading in the measured roughness 
exponents, scaling anisotropy were looked for [63, 62, 64] in the fracture surfaces of 
various materials (quasi-crystals, glasses, metallic alloys, wood and mortar) broken 
under various conditions (cleavage, stress corrosion, tension) and various velocities, 
ranging from 10~ 12 m/s to 10 2 m/s. At this point, it is important to set the 
notations (Figure 14): In all the following, the axis e x , e y and e z will refer to the 
direction of crack propagation, the direction of tensile loading and the mean direction 
of the crack front, respectively. Special attention was paid to ensure that fracture 
surfaces were observed in a steady state regime, far enough from crack initiation so 
that roughness becomes statistically invariant of the position x along the direction 
of crack growth. By computing the height-height correlation function Ahg(Ar) for 
profiles making various angles 9 with e x , the existence of anisotropy in the scaling 
features (Figure 7) was evidenced: £(#) varies from = 0) = C|| — 0.6 (profiles taken 
parallel to the direction of crack growth) to ((6 = 90°) = C-L — 0.8 (profiles taken 
perpendicular to the direction of crack growth, parallel to the mean direction of crack 
front). 

As a consequence, a full characterization of the scaling statistical properties of 
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Figure 7. (a) ID height-height correlation function measured parallel and 
perpendicular to the crack propagation direction. The axes are logarithmic, (b) 
Variation of roughness exponent £ as a function of the angle 8 between the 
considered profile and the direction of crack propagation. Data were measured 
in quasi-crystals (From [62]). 



fracture surfaces calls for the computation of the 2D height-height correlation function: 

Ah(Az, Ax) = ((h(z + Az,x + Ax) - h(z, x)f) z , x (16) 

This was shown [63, 62, 64] to take a peculiar form, referred to as Family- Viseck 
scaling [ ]: 



Ah(Az, Ax) oc Ax^f(Az/Ax^ /c± ) 
where /(it) oc \ V, •!• 

J K ' y vS 1 - if u» 



(17) 

if u< c v ; 



c 

Such Family- Vicsek scalings are classically observed in interface growth processes. 
In this respect, the exponent C|| an d the ratio Ci/C|| can be mapped to a growth 
exponent and a dynamic exponent, respectively. It is worth to mention that such 
an anisotropic Family- Viseck scaling of fracture surfaces was first suggested by the 
phcnomcnological Langevin description of crack propagation proposed by Bouchaud 
et al (1993) [93] detailed thereafter in section 6.3. 

The two exponents (± — 0.75±0.05 and C|| = 0.6±0.05 were found to be universal, 
independent of the considered material, the failure mode and the crack growth velocity 
over the whole range from ultra-slow stress corrosion fracture (picometer per second) 
to rapid failure (few hundreds meter per second). On the other hand, the range of 
length-scales over which these Family- Viseck scaling features are observed is limited 
and material dependent. This will be discussed in section 4.3. 

By rescaling the distances Aa; and Az by the topothesies l x and l z defined as 
the scales at which the out-of-planc increment becomes equal to the in-plane one, i.e. 
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Figure 8. Family- Viseck scaling of the 2D height-height correlation function in 
quasi-crystals (From [62]) and in mortar (From [ ]) using equation 17. 



Ah(Az = 0,Ax = £ x ) — £ x and Ah(Az = £ z ,Ax = 0) = £ z , one can rewrite the 
Family- Viseck scaling [62]: 



Ah(Az, Ax) =£ x (Ax/£ x ) c i\g 



Az/f, 



where g(u) 



(18) 
1 if u< 1 
u c± if u> 1 



The form of / is then found to be universal [6 ] , independent of the considered material, 
of the failure mode and of the crack growth velocity (Figure 9) . 

It is important to insist that this universal Family- Viseck scaling is observed 
in fracture surfaces far enough from crack initiation, in a regime where roughness 
becomes statistically invariant with respect to translations along the direction of crack 
growth e x . The transient roughening regime starting from an initial straight notch, 
was extensively studied by Lopez, Morel et al [91, 75, 94]. It was found to display a 
scaling significantly more complex than that expected in Family- Viseck scenario. It 
involves the global roughness exponent C g i b defined at the end of section 4.1 as well 
as a new independent dynamic exponent. These two exponents are found to depend 
on both the material and specimen geometry. 



4-.S. On the relevant length- scales 

Several fractography observations reported in the literature are found not to be 
compatible with the Family- Viseck scaling (equation 17) and its associated exponents 
— 0.75, £|| ~ 0.6}. In particular, experiments in sandstone [77, 78, 66], in artificial 
rocks [S ] and in granular packings made of sintered glass beads [80, 18, 81] have 
shown self-affine scaling properties characterized by roughness exponents significantly 
smaller, around 0.4—0.5. In the case of granular packings made of sintered glass beads, 
the 2D height-height correlation function was found to exhibit Family- Viseck scaling 
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Figure 9. Collapse of the dimensionless 2D height-height correlation function 
using Equation 18 in various materials, (From [62]). 
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Figure 10. Variation of the 2D height-height correlation function Ah with Az for 
various values of Ax for granular packings made of sintered glass beads. The inset 
shows the direct curves while the data collapse in the main graph was obtained 
from equation 17 using = 0.4 and £| = 0.5 (From [18]). 



as for the materials investigated in the previous section, but with {(± ~ 0.4, 0i ~ 0.5} 
[18]. This suggests the existence of a second universality class for post-mortem fracture 
surfaces. 

To uncover the origin of these two distinct series of observed exponents, The range 



CONTENTS 



18 



of length-scales over which these two regimes are observed was examined [18]. In all 
experiments compatible with {(± ~ 0.8, £|| ~ 0.6}, fracture surfaces were observed at 
small length-scales, below a given cutoff length £ of e.g. few nanometers in quasicrystals 
[62], few tens of nanometers in glass [73, 18, 95], few hundreds of micrometers in 
metallic alloys [73, 72], few millimeters in wood [75, 63], few tens of millimeters in 
mortar [79, 94]... On the other hand, observations compatible with {(± ~ 0.4, £|| ~ 0.5} 
were observed at large scales, well above the microstructure scale, up to a cutoff length 
set by the specimen size [66]. 

It was then proposed, in ref. [18], that the upper cutoff length £ that limits the 
{(± ~ 0.8, 1 ~ 0.6} Family- Viseck scaling regime is set by the size of FPZ. This 
conjecture was proven to be true in oxide glasses [18], quasicrystals [62] and mortar 
[94]. Work in progress [96] seems to indicate that it is also true in metallic alloys: CT 
specimen were broken at various temperatures, ranging from 20 K to 148 K - this 
allows to tune the FPZ size, from 20 /im to 1 mm. And in all these tests, the cutoff 
length £ of the self-affine £ ~ 0.8 regime is found to be roughly proportional to the 
process zone size. 

It is worth to mention that no anomalous scaling was reported in Sandstone 
fractography experiments [ ]. This suggests that the anomalous scaling commonly 
associated with the £ ~ 0.8 self-affine regime [91, 75] is absent in the large scale 
£ ~ 0.4 self-affine regime. It is then appealing to associate the observation of this 
anomalous scaling to damage spreading at scale smaller than the FPZ. 

Let us note also that a third isotropic scaling regime arising at small scale, 
characterized by a roughness exponent £ ~ 0.5 was observed in ductile materials like 
Ti 3 Al-based super-a2 intermetallic alloy broken under fatigue within 1 nm — 10 /im 
range [73, 10, 97], and in Zr-based metallic glass within the range 100 nm — 1 mm [98]. 
It should be emphasized that in the second case, the distribution of height increments 
does not display the collapse given by equation 13 expected for truly self-affine profiles. 
Similar £ ~ 0.5 isotropic regime was also claimed [73, 10] to be observed on the AFM 
fracture surfaces of soda-lime broken under stress corrosion, within 1 nm— 10 nm range. 
However, this observation has been recently called into question since the observations 
remain confined within (in-plane) AFM resolutions. Moreover, it was not reproduced in 
nanoresolved fracture surfaces of various silicate glasses (silica, soda-lime, borosilicate 
and aluminosilicate) broken under stress corrosion with crack velocities as small as the 
picometer per second [63, 18]. As a result, this third small-scale regime characterized 
by £ ~ 0.5 is suspected to be inherent to ductile fracture. As such, it will not be 
discussed in the sequel. 

Let us conclude this section by drawing the outlines of the fractography 
experiments reviewed here: 

(i) Failure of brittle materials leads to rough fracture surfaces, with self-affine 
morphological features. 

(ii) Far enough from crack initiation, the morphology of fracture surfaces exhibits 
anisotropic Family- Viseck scaling, characterized by two distinct roughness 
exponents £|| and £_l, measured parallel and perpendicular to the direction of 
crack propagation, respectively. By analogy with interface growth problems, the 
exponent £|| and the ratio Cl/C|| can be mapped to a growth and a dynamic 
exponent, respectively. 

(iii) Two distinct sets of universal exponents are observed: {£j_ ~ 0.8, £|| = 0.6} at 
small scales and {Cj_ — 0.4, £|| = 0.5} at large scale. 
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Figure 11. (a) Random Fuse Networks with a square lattice geometry. The 
system is "loaded" either by imposing the voltage U at the upper electrode 
(fracture test with imposed displacement conditions) or by imposing the total 
current / crossing the system (fracture test with imposed force conditions), (b) 
Each fuse has a unity conductance and a randomly assigned breakdown threshold 

{Uc,Ic}. 



(iv) The cutoff £ of the small-scale regime is set by the FPZ size. 
5. Numerical observations in lattice models 

The two preceding sections allowed to illustrate that the failure of brittle heterogeneous 
materials leads to (i) intermittent crackling dynamics and (ii) rough fracture surfaces. 
In both cases, one observes scaling invariance in term of distribution and morphology. 
Some of these scaling features appear universal - the associated scaling exponents are 
the same in a wide range of materials and loading conditions. This is reminiscent of 
critical phenomena where the balance between local (thermal or quenched) disorder 
and simple interactions between microscopic objects may select such complex scale 
invariant organization of the system (see e.g. [99] for an introduction to critical 
phenomena) . 

In this context, it was proposed, initially by de Arcangelis et al (1989) [100], 
to exploit the analogy between scalar mode III elasticity and electricity and sketch 
an heterogeneous brittle material as a network of fuses of unit ohmic resistance and 
randomly distributed breakdowns (Figure 11). In these so-called Random Fuse Models 
(RFM), current and voltage are analogous to force and displacement, respectively. The 
goal is not so to reproduce exactly the failure of real brittle materials, but to see to 
which extent one can reproduce the scaling features presented in Sections 3 and 4 
in such minimalist materials, keeping only the two main ingredients characterizing 
material failure, namely the microstructure randomness and the long-range coupling 
that accompanies the load redistribution after each element breakdown. 
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Figure 12. Successive snapshots in a typical simulation of 2D 100 X 100 RFM 
(square lattice) "loaded" by gradually increasing the voltage U. The red large 
points correspond to the cascade of fuses breakdown at the considered value of U . 
The gray small points corresponds to the fuses that have been burnt from (7 = 
to the considered value of U . In the last image of this sequence, the lattice is 
broken into two pieces. Note the roughness of the fracture line. 



5.1. Intermittency 

A typical "fracture" experiment performed on 2D RFM is represented in Figure 12. 
The "loading" is increased by increasing the voltage U of the nodes belonging to the 
upper electrode while keeping to zero the bottom ones. Voltages at each node and 
currents crossing each element are determined such that Kirchoff law and Ohm's law 
are satisfied everywhere. At some point, a fuse breaks. As a result, the current increases 
abruptly in the remaining elements. In turn, this may trigger a cascade, avalanche, 
of elements breakdowns that can either lead to a stable situation or to the failure of 
the whole structure. Typical evolution of the size of these avalanches as a function 
of the imposed voltage U is plotted in figure 13. Avalanches are analogue to the AE 
events observed prior to the catastrophic failure of real brittle materials, as discussed 
in section 3.2. 

The distribution P(S) of avalanche sizes, i.e. the number of fuses participating in 
a breakdown cascade has been investigated numerically, first by Hansen and Hemmcr 
(1994) [101] and later by other teams [102, 103, 104, 105]. It was found to follow a 
power-law P(S) oc S~ T up to a cut-off So that scales with the lattice size L as So oc L D 
where D refers to the avalanche fractal dimension. Initially, 2D simulations yielded 
a universal value r ~ 2.5 [101, 102]. However, recent large-scale simulations suggest 
that r depends slightly on the lattice type in 2D systems: r = 2.75 and t = 3.05 in 
diamond and triangular lattice respectively [104]. On the other hand, D was found to 
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Figure 13. Evolution of the avalanche size S, i.e. the number of fuse participating 
in a breakdown cascade as a function of the applied voltage U in the RFM 
simulation presented in Figure 12. 



be universal: D ~ 1.18. In 3D simulations, both exponents were found to be universal 
[105]: t = 2.5 and D = 1.5. 

5.2. Morphology of fracture surfaces 

The morphology of final cracks has also been investigated numerically [106, 107, 108, 
103, 109, 104, 110, 111]. As in experiments, it was found to exhibit self-affine scaling 
features. In 2D RFM such as the one shown in Figure 12, large scale simulations 
yield a universal roughness exponent £ = 0.71 independent of the lattice topology, the 
presence or not of a pre-existing initial notch and of the disorder in breaking thresholds 
[104, 111]. In 3D systems, C - 0.42 [108, 103, 110]. 

It is interesting to note that the value of £ measured in 2D RFM is pretty close to 
the C — 0.7 observed experimentally in quasi-two-dimensional materials [83, 84, 65]. 
In 3D RFM, the value £ ~ 0.42 are significantly smaller than the £ ~ 0.8 observed 
experimentally in a wide range of materials (see sections 4.1 and 4.2). On the other 
hand, it is pretty close to the value £ ~ 0.4 — 0.5 observed at scale larger than 
the FPZ size in sandstone [77, 66], Glassy ceramics made of sintered glass beads 
[80, 81] and mortar [94]. It should be emphasized that the morphological scaling 
analyses reported in [108, 103, 110] were performed by averaging together the results of 
different simulations starting from configurations without initial notch, and therefore 
no prescribed direction of crack propagation. As a consequence, the scaling anisotropy 
reported in experiments and described in section 4.2 cannot be captured. 

Let us finally mention that large scale numerical simulations revealed the existence 
of anomalous scaling in both 2D and 3D RFM [104, 111]. The origin of this anomalous 
scaling remains unclear. It seems to be intrinsic to the scalar dimensionality of RFM 
since it disappears in Random Beam Models [112]. 
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Figure 14. Sketch and notation of a crack front propagating in a 3D material. 



6. Stochastic theory of crack growth 

Experimental observations gathered in sections 3 and 4 revealed that slow crack growth 
in brittle disordered materials (i) exhibits an intermittent dynamics with scale free 
distributions in space, time and energy and (ii) leads to rough fracture surfaces with 
scale invariant morphological features. As was exposed in Section 5, these two aspects 
can be qualitatively reproduced in simple numerical models such as RFM. However, 
these models rely on important simplifications which make quantitative comparisons 
with experiments difficult. 

We turn now to another approach - pioneered by Gao and Rice (1989) [ ] and 
later extended by Bower and Ortiz (1991) [113], Schmittbuhl et al (1995) [16], Larralde 
and Ball (1995) [114], Ramanathan et al (1997) [115] and Lazarus (2003) [116] which 
consists in using LEFM to describe crack growth in an elastic isotropic material in 
presence of deterministic or random obstacles. As we will see below, to first order, 
the equation of motion and the equation of path are decoupled and can be solved 
independently. 

6.1. Equation of motion 

Let us consider the situation depicted in Figure 14 of a crack front that propagates 
within a 3D elastic solid. Provided that the motion is slow enough, the local velocity 
of a point M(z,x = f(z,t),y — h(z,x — f(z,t))) is proportional to G(M) — r(M) 
(Equation 5). Obstacles and inhomogeneities in the material structure are then 
captured by introducing a fluctuating component into the fracture energy: T(M) = 
r° (1 + rj(M.)) where r)(M.) denotes the fluctuating part of fracture energy. This 
induces distortions of the front, which in turn generates perturbation in G(M). One 
can then use the three-dimensional weight-function theory derived by Rice[117] to 
relate G(M) to {/, h}. To first order, G(M) depends on the in-plane distortions f(z, t) 



CONTENTS 



23 



only and the resulting equation of motion reads [16, 115]: 




dz' + rVz,f(z,t)), 



(19) 



where F = G° — T° and G° denotes the reference mechanical energy release which 
would result from the same loading with a straight front at the same mean position. 
As it, G° depends on the macroscopic geometry and loading conditions. It is worth to 
mention that this equation of motion was recently shown by Dalmas and co-workers 
[118] to reproduce quantitatively experiments which consists in making a planar crack 
propagate along a weak interface with simple mono-dimensional regular patterns in a 
thin film stack (single band or periodic bands drawn parallel to the direction of crack 
propagation) . 

In the case of heterogeneous materials, ?y(M) is chosen to be random, with short 
range spatial correlation, zero mean and constant variance. Let us first consider 
the situation of constant remote stress loading conditions, i.e. constant F - This 
assumption will be released later in this section. As first noticed by Schmittbuhl et al 
(1995) [16], the front motion described by Equation 19 exhibits a so-called depinning 
transition controlled by the "force" F = G° — T° and its position with respect to a 
critical value F c (see Figure 15a). 

• When F < F c the crack front is pinned by disorder and does not propagate. 

• When F 3> F c the crack front moves with a mean velocity v = (df/dt) Zyt 
proportional to F: v / fi = G a — T°. One recovers the equation of motion expected 
by LEFM (Equation 5). 

• When F = F c a critical state is reached. 

The consequences of this criticality will be extensively commented later in this 
section. But it is now time to discuss more explicitly the form taken by G° and its role 
in the intermittent crackling dynamics discussed in Section 3. In this respect, it has 
been suggested, in ref. [17], to consider the case of a crack growing stably in a material 
remotely loaded by imposing a constant displacement rate. This situation is the one 
encountered as e.g. in earthquakes problems where a fault is loaded because of the 
slow continental drift, or in the interfacial experiments described in Section 3.3 where 
a crack front is made propagate along the weak heterogeneous interface between two 
Plexiglas block by lowering the bottom part at constant velocity. In this situation, G° 
is not constant anymore, but: 

• increases with time since the remote loading increases with time. 

• decreases with the mean crack length (f)(t) = (f(z,t)) Zyt since the material 
compliance decreases with (/) (t) . 

As a result, provided that the mean crack growth velocity is slow enough and the 
mean crack front is large enough, one can write [17]: 



where c and k are constants depending on the precise geometry and the loading 
conditions (see e.g. [17] for their value in the case of the interfacial experiment 
described in Section 3.3). The crack motion can then be decomposed as follows (Figure 
15b): When F < F c , the front remains pinned and the effective force increases with 
time. As soon as F > F c , the front propagates, (f)(t) increases, and, as a consequence, 



F(t, (f)) = ct-k(f)(t), 



(20) 
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Figure 15. (a) Pinning/depinning transition predicted by Equation 19. The 
control parameter in this system is the force F defined as the difference between 
the stress intensity factor K® that would result from the same loading with 
a straight crack front at the same mean position and the mean value of the 
local toughness K®. The order parameter is the mean value of the velocity 
v = (df/dx) z t averaged over all z and t. (b) Time evolution of F(t) in 
displacement imposed loading conditions (Obtained by solving numerically Eq. 
19 with fj, = G° = r° = 1 and F(t, (/>) given by Eq. 20 with c = 10 -5 and 



F is reduced. This retro-action process keeps the system attracted to the critical state, 
as in models of self-organized criticality [119]. For high values of c, the front moves 
smoothly with mean velocity (v) — c/k. For for low values of c, stick-slip motion is 
observed and the crack propagates through distinct avalanches between two successive 
pinned configurations. In the limit c — > and k — > 0, criticality is reached and one 
expects to observe universal scaling features: 

• independent of the microscopic and macroscopic details of the system, i.e. the 
considered material and the loading conditions 

• identical to those observed in other systems belonging to the same universality 
class, as e.g. interface motion in disordered magnets [120, 121] and wetting of 
rough substrates [122, 123]. 

These scaling features can be predicted using theoretical tools issued from 
Statistical Physics like Functional Renormalization Methods (FRG) or numerical 
simulations. In particular, for F slightly above F C1 the front f(z) is self-affine and 
exhibits Family- Viseck dynamic scaling form up to a correlation length £: 



where Qh and k refer to the roughness exponent and the dynamic exponent 
respectively. It is worth to recall that £h is associated to the in-plane roughness of 
the crack front and therefore is different from £ defined in section 4 to characterize 
the out-of-plane roughness. Above £, the front is no longer self-affine and exhibits 
logarithmic correlation. 



k = 10~ 3 ). 



((f(z + Az,t + At) -f(z,t)) 2 )^ 2 cx At^/ K g(Az/At^ K ) 




(21) 



where 
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The exponents and n were estimated using FRG methods. They were found 
to be between {( H = 1/3, k = 7/9} [ ] (first order) and {( H = 0.47, k = 0.66} [124] 
(second order). More recently, they were precisely evaluated by Rosso and Krauth 
[125] and Duemmer and Krauth [126] using numerical simulations: 

( H = 0.385 ±0.005, k = 0.770 ±0.005 (22) 

The correlation £ is set by the distance between the force F and its critical value F c . 
In situations where F is kept constant, £ diverges at F c as: 

£ oc (F- F c )~ v with v = - 1 = 1.625 ±0.010, (23) 

i — Ch 

where the relation between i/ and comes from a specific symmetry of Equation 19, 
referred to as the statistical tilt symmetry [127]. In situations where F varies with time 
and crack length according to equation 20, £ is set by the "strength" of the feedback 
term. To be more precise, £ is defined [121] as the length for which the feedback term 
-k(f)(t) in Equation 20 balances the elastic term (G°/2ir) J (f(z)-f{z'))/{z-z') 2 dz' 
in Equation 19. This leads to kf^/L w G //£ where L refers to the specimen width. 
Hence, equation 23 should be replaced by [121]: 

foe AT"* with i/ fc = 1/2 (24) 

The front propagation occurs through avalanches between two successive pinned 
configurations (Figure 16a). An avalanche of size S translates into an increment S/L 
for the mean crack length (Figure 16b). Equation 20 ensures also that the mechanical 
energy E released during an avalanche is also proportional to S (Figure 16c). As a 
result, the energy signal is very similar to that observed in figure 2a. 

The avalanche size S is then shown [123] to obey power-law distribution: 

P(S) oc S- T f(S/S ) with t = 2 - 1 + 2l/ f g = 1.280 ±0.010, (25) 

K 1 + Ch) 

where the cut-off scales as So oc £ 1+ ^ H , i.e., after having used Equation 24: 

S oc k- 1/a with g = — 2 — - = 1.445 ± 0.005 (26) 

1 + Ch 

The avalanche duration T scales with the avalanche size I as T oc £ K . Since the 
avalanche size S scales as i 1+< * H , one gets: 

T oc S a with a = — ^— = 0.555 ± 0.005 (27) 
1 + Ch 

And, after having changed the variables accordingly in equation 25, one shows that T 
is distributed as: 

P{T) oc T~ a f(T/T ) with a = 1 + V -— - = 1.500 ± 0.010, (28) 

KV 

where the cut-off scales as To oc £ K , which implies: 

T oc fc" 1/A with A = l/nv k = 1.300 ± 0.010 (29) 

Presently, these scaling predictions on the global dynamics of the crack front were 
not confronted directly to experiments. Work in progress [128] seems to indicate that 
the size and duration of the mean front jumps observed in the interfacial experiment 
presented in section 3.3 are distributed according to equations 25 and 29, respectively. 
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(a) 




Figure 16. (a) Sketch of the front propagation between two pinned configuration. 
The avalanche area A scales with the size I of the depinned zone as A oc i 1Jr ^ H . 
(b) Typical evolution of the mean crack length. The crack progresses through 
jumps of size S proportional to A. (c) Typical evolution of the energy released 
E. The energy released during an avalanche is also proportional to S. Note the 
similarity between this curve and the one presented in figure 2a 



On the other hand, the local dynamics of the stochastic description presented here 
was directly confronted to the intcrfacial experiment presented in section 3.3 [18]. The 
equation 19 with F given by equation 20 was solved for various values of parameters 
and the results were analyzed using the procedure described in section 3.3: The map 
of waiting times (figure 17a), i.e. the time spent by the crack front in each point (z, a;) 
of the recorded region was computed, as in the experiments presented in section 4.3. 
Then, the avalanches were defined by thresholding these maps of waiting time, also 
as in section 4.3 (figure 17b). The statistics of the obtained avalanches was finally 
computed and found to be the same as that observed in the interfacial experiment, 
with in particular a power-law distribution for avalanches area characterized by an 
exponent t w ~ 1.7 (figure 17c). Presently, the relation between r w and the standard 
exponents defined in equations 22 to 29 remains unknown. 

It is finally worth to discuss the roughness exponent (h — 0.385 expected for the 
in-plane crack front (equation 22). This value is significantly smaller than the value 
~ 0.6 reported in interfacial experiments [86, 87]. However, recent experimental 
results [129] shed light on the apparent disagreement between the two and report at 
large scales a roughness exponent (h — 0.35 compatible with theoretical predictions 
of equation 22. 
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Figure 17. (a) Typical gray scale map of the waiting time matrix w(x, y) obtained 
from the solution of equations 19 and 20. (b) Spatial distribution of clusters 
corresponding to velocities four times larger than the mean one. The clusters 
duration is given by the clusters color according to the colorscale given in inset. 
The distribution of the cluster size S are plotted in (c). The various symbols 
corresponds to various values of c, k and cluster thresholds. The straight lines 
correspond to P(S W ) oc S w Tm with t w = 1.7 (From [17]). 



6.2. Path equation 

We turn now to out-of-plane crack roughness. The universal morphological scaling 
features observed both in experiments and in lattice simulations reported in sections 
3 and 4.2, respectively, have lead to intense theoretical development over the last 
20 years. We will focus here on the stochastic descriptions of crack path in 3D 
elastic disordered materials, in the spirit of what was done in the preceding section 
about motion equation, returning in section 6.3 to a brief discussion of alternative 
interpretations of the morphological scaling features. 

Let us consider again the situation depicted in figure 14 where a crack front 
propagates within a 3D isotropic elastic solid remotely loaded in mode I. Provided 
that the motion is slow enough, PLS imposes that the path chosen by the crack is 
the one for which the net mode II stress intensity factor vanishes at each point M of 
the front at each time t (section 2.3): Ku(M.(t)) = 0. The tough part of the problem 
is to relate Kjj(M.(t)) to the crack roughness f(z,t) and h(z,x) (figure 14). This 
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has been made possible thanks to the work of Larralde and Ball (1995) [ ] later 
refined by Movchan et al (1998) [130] which shows that, to first order, Kjj(M.(t)) 
only depends on h. It can then be written as the sum of two contributions: A first 
contribution arises from the coupling of the singular mode I component Kj of the 
stress field of the unperturbed crack with the position along the crack edge and a 
second contribution comes from the coupling between the slope of the crack surface 
and the non singular T° normal stress in the direction of crack propagation. Both 
contributions can be computed analytically (see [ ] for the complete form) , but the 
second contribution was shown to be negligible compared to the first one when one 
considers roughness length-scales small with respect to the system size [114]. To the 
remaining contribution, a stochastic 3D term is added to capture the heterogeneous 
nature of the material, and a constant term Kjj is added to take into account the 
inherent misalignment in any loading system. Introducing all these terms in the PLS 
Kjj(M.(t)) = leads to the following path equation [18]: 



At this point, it is worth to recall that this path equation was derived within 
elastostatic and, as such, cannot work for dynamic fracture when the crack speed 
exceeds a few percent of the Rayleigh wave speed. 

The form taken by this equation is very similar to that of Eq. 19 which describes 
the in-plane motion of the crack front. Two main differences are however worth to be 
discussed: 

(i) The absence of explicit time t in equation 30. This comes from the fact that to 
first order, Kjj only depends on the out-of-plane roughness h(z,x). This yields 
a decoupling between this path equation and equation 19 describing the crack 
dynamics. It implies that the scaling properties of fracture surfaces will not depend 
on the crack velocity v or crack loading G°, contrary to those of the in-plane 
projection of the crack front. 

(ii) The 3D nature of the stochastic term rj(z^ x, h(z, x)). 

This second point makes equation 30 extremely difficult - if not impossible - to solve. 
Two limit cases can however be considered. 

The first limit consists [114, 115] in considering very smooth fracture surfaces so 
that the 3D stochastic term r)(z,x,y = h(z,x)) reduces to an effective "thermal" 
term rj(z,x). Fracture surfaces are then predicted to exhibit logarithmic scaling. 
This prediction seems in apparent disagreement with the experimental observations 
reported in section 4. However, recent observations performed by Dalmas and 
coworkers (2008) [ ] on nanoscale phase separated glasses reveals logarithmic 
roughness at large length scales. 

The second approximation consists [18] in writing ij(z,x,y — h(z,x)) as the sum 
of two 2D terms: ri{z,x,y = h(z,x)) = ri(z,x) + rj(z,y = h(z,x)). The morphology 
of the fracture surface h(x, z) is then given by the motion of the elastic string 
h(z) that "creeps" - the x coordinate playing the role of time - within a random 
potential r/(z,h(z,x)) due to the "thermal" fluctuations i](z,x). The fracture surface 
h(z,x) obeys then to Family- Viseck anisotropic self-affine scaling (equation 17) with 

— 0.385, C|| — 0.5} in perfect agreement with observations reported at large scales, 
as e.g. in sandstone [77, 78, 66], in artificial rocks [ ] and granular packing made of 
sintered glass beads [80, 18, 81] (see also section 4.3). 




(30) 
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Let us finally mention that similar stochastic approach has been developed 
recently by Katsav et al [131] in simpler 2D elastic materials. In this case, the path 
equation can be solved analytically and leads to apparent self-affine scaling features 
at small scales, with £ = 0.5, and departure from self-affinity at large scales. However, 
these predictions seem in contradiction with experimental observations [83, 84, 65] 
which report a unique £ ~ 0.7 self-affine scaling regime in quasi- two-dimensional 
materials. 

6.3. Alternative interpretations for self-affine crack roughness 

The stochastic theory of crack growth in elastic brittle disordered materials described 
in section 6.2 seems to reproduce quite well the morphological scaling properties 
of cracks for length-scale above the FPZ size. On the other hand, it cannot 
describe observations performed below the FPZ size. In particular, in essence, it 
cannot reproduce the small-scale anisotropic Family- Viseck scaling characterized by 
a roughness exponent £ ~ 0.8 reported in a wide range of materials and discussed in 
section 4. In this respect, we will present now some alternative approaches that were 
developed to account for the morphological scaling features observed in fractography 
experiments as well as in RFM simulations. 

Hansen et al (1991) [106] have studied the morphology of crack lines (2D case) 
or crack surfaces (3D) obtained after the breakdown of RFM with perfectly plastic 
fuses, i.e. fuses that act as unit Ohmic resistors up to a threshold voltage and that 
carry a constant current above. They have shown that the problem can be mapped 
to that of a directed polymer (2D case) or a minimum energy surface (3D case) in 
random potential. As a result, crack path can be mapped to a KPZ equation [ ] and 
crack lines/surfaces are expected to display self-affine morphological isotropic scaling 
features with a universal roughness exponent equal to £ = 2/3 in 2D, and £ ~ 0.42 in 
3D. They conjectured that these results, demonstrated rigorously on elastic-perfect- 
plastic model, can also be relevant for brittle fracture. Indeed, these predictions are in 
close (but not perfect) agreement with the £ ~ 0.71 (2D) and £ ~ 0.42 (3D) obtained 
from numerical simulations of brittle RFM and presented in section 4.2. However, the 
anomalous scaling is not recovered. Hansen et al 's predictions are also very close 
to the £ ~ 0.7 observed experimentally in 2D sheet of paper [83, 84, 65] and the 
£ ~ 0.4 — 0.5 observed in sandstone [77, 78, 66] or glassy ceramics made of sintered 
glass beads [80, 18, 81]. 

Bouchaud et al (1993) [93] proposed to model the fracture surface as the trace 
left by a line moving through a 3D disordered landscape - the dynamics of which is 
described through a phenomenological nonlinear Langevin equation derived initially 
by Ertas and Kardar [133] to describe the motion of vortex lines in superconductors. 
The resulting fracture surfaces are then predicted to display anisotropic Family- Viseck 
scaling with {£j_ ~ 0.75, £h ~ 0.5}. This first work was then extended by Daguier et 
al (1997) [73] who showed that, depending on the crack velocity v, two distinct self- 
affine regimes can be expected: A large length-scale regime with {£j^ ~ 0.75, £n ~ 0.5}, 
and a small length-scale regime with {£^ ~ 0.5, £m ~ 0.4}. The crossover £ between 

the two is expected [73, 10] to diverge with v as £ oc The exponents 

values obtained at large scale/large velocities are in excellent agreement with those 
observed in many brittle materials at scales smaller than the FPZ size (see section 
4.3). However, the existence of a second small scale regime all the more important 
than crack velocity is large is incompatible with experiments. 
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More recently, Hansen, Schmittbuhl and Batrouni (2003) [134, 135] suggested 
that the universal scaling properties of fracture surfaces are due to the fracture 
propagation being a damage coalescence process described by a stress-weighted 
percolation phenomenon in a self-generated quadratic damage gradient. 

Then, assuming that the damage profile is proportional to the system size L, the 
rms roughness w of the crack line /surface is expected to scale as: 

maL^ with £ = — , (31) 

1 2z^p erc 

where v perc is the exponent describing the divergence of the correlation length 
£, P erc <x (p ~ Pc) ~ VpaTC as the density of broken bonds p is reaching the critical value 
p c in standard percolation problems (see e.g. [136] for an introduction to percolation 
theory) . 

In percolation theory, the exponent v perc depends only on the system dimension 
and dimensionality. In standard scalar percolation, v perc = 4/3 in 2D, and v perc — 0.88 
in 3D [136]. This leads to ( = 8/11 in 2D, and C ^ 0.64 in 3D. The first value is 
compatible to the £ ~ 0.7 measured in 2D brittle RFM (see section 4.2 or references 
[104, 111]), while the second one is significantly higher than the £ ~ 0.42 observed 
in 3D brittle RFM (see section 4.2 or references [108, 103, 110]). In tensorial elastic 
percolation, v perc = 2 in 3D, which leads to £ = 4/5. The value of the roughness 
exponent is in excellent agreement with £ ~ 0.8 observed in many brittle materials 
observed at scales smaller than the FPZ size (see section 3), but the Family- Viseck 
anisotropic scaling is not captured. 

This interpretation yields recent controversies. First, the analysis of large-scale 
RFM simulations [137] suggests that, contrary to what was assumed, (i) damage 
profile is not quadratic and (ii) the width of the damage profile is not linear with 
the system size. Second, it was mentioned [138] that the exponent £ in equation 31 
cannnot be intepreted as a "standard" roughness exponent since, strickly speaking, 
the front is not self-affine. Indeed, in such a model, the fracture perimeter displays 
substantial overhangs the size of which is comparable with their width. To be more 
precise, self-similarity (i.e. isotropic scaling) is expected for length scales smaller than 
to, and a trivially flat front beyond. In response, Hansen, Schmittbuhl and Batrouni 
[139] pointed out that experimental fracture fronts are also found to display significant 
overhangs while remaining self-affine at large scales. They argue that the non-isotropic 
scaling given by equation 31 is sufficient to define self-affine surfaces. 
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7. Concluding discussion 

Many experiments and fields observations have revealed that brittle failure in materials 
exhibit scale-invariant features (see Tab. 1 for a summary). In particular: (i) Fracturing 
systems displays jerky crackling dynamics with random impulsive energy release, 
as suggested from the AE accompanying the failure of various materials and, 
at much larger scale, the seismic activity associated to earthquakes. The energy 
distribution of these discrete events forms a power-law that spans over many orders of 
magnitude, with no characteristic scale and (ii) roughness of cracks exhibits self-affinc 
morphological features, characterized by roughness exponents. These observations are 
common to many brittle materials and can be reproduced qualitatively in the electrical 
breakdown of Random Fuse Network. On the other hand, they cannot be captured by 
standard LEFM continuum theory. 

Some of these scale-free distributions and scale-invariant morphological features 
are universal. Others are not. The experiments and lattice simulations presented in 
sections 3, 4 and 5 allows to distinguish two cases: 

(A) Damage spreading processes within a brittle material preceding the initiation of a 
macroscopic crack. The associated micro-fracturing events release energy impulses 
which are power-law distributed. The associated exponent is non-universal, but 
depends on the considered materials, loading conditions environment parameters 
(see section 3.2)... This transient damage spreading is also suggested [64, 94] 
to be responsible for the anomalous non-universal scaling exhibited in the 
initial transient roughening regime of fracture surfaces following crack initiation 
[91,75,94]. 

(B) Macroscopic crack growth within a brittle material. When this propagation is slow 
enough, it exhibits an intermittent crackling dynamics, with sudden jumps and 
energy release events the distributions of which form power-law with apparent 
universal exponents. Crack growth leads to rough fracture surfaces that display 
Family- Viseck universal scaling far enough from crack initiation. 

Quite surprisingly, seismicity associated with earthquakes seems to belong to the 
second case and exhibits quantitatively the same statistical scaling features as that 
observed in experiments of interfacial crack growth along weak disordered interfaces 
[58]. 

The (non-universal) scaling features associated with the microfracturing events 
observed experimentally in case (A) are reproduced qualitatively in lattice models such 
as the RFM presented in section 4. Nevertheless, what sets precisely the distribution 
and the value of scaling exponents in experiments remains largely unknown. On the 
other hand, the recent stochastic extensions of LEFM theory presented in section 6 
seems efficient to capture quantitatively the observations performed in case (B). In 
these stochastic LEFM descriptions, the onset of crack growth is found to be analogous 
to a critical transition between a stable phase where the crack remains pinned by the 
material heterogeneities and a moving phase where the mechanical energy available 
at the crack tip is sufficient to make the front propagate [16, 115]. While growing, 
the crack decreases its mechanical energy and gets pinned again. Provided that the 
growth is slow enough, this retro-action process keeps the system close to the critical 
point during the whole propagation, as for self-organized-critical systems [140]. As 
seen in section 6.1, the resulting dynamics exhibits spatio-temporal intermittency 
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characterized by universal statistical features that can be predicted theoretically using 
FRG methods and are compatible with experimental observations. 

The analogy between the onset of crack propagation and a critical dynamic 
transition has numbers of other consequences. Roux, Charles et al (2003,2004) 
[141, 142] made use of the universality manifested around the depinning onset 
to determine the distribution of effective macroscopic toughness in heterogeneous 
materials. Their main results are: (i) Universal scaling between toughness variance 
and specimen size, and (ii) the existence of a universal toughness distribution for 
large enough specimens. These predictions were shown [143] to reproduce fairly well 
the statistics of crack arrest lengths observed in indentation experiments performed in 
various brittle materials. More recently, Ponson et al (2007,2009) [144, 145] investigate 
the role of a finite temperature in this stochastic LEFM description of crack growth 
and propose a creep law to relate crack velocity and stress intensity factor in sub- 
critical failure regime. This was shown recently to describes rather well experiments 
of paper peeling [49] and subcritical crack growth in sandstone [145]. 

It should be emphasized that the stochastic descriptions presented in section 
6 were developed within elastostatic approximation. As such, they cannot account 
for the dynamic stress transfers through acoustic waves occurring as a dynamically 
growing crack is interacting with the material disorder [146, 147, 148, 149, 150]. The 
recent availability of analytical solutions for weakly distorted cracks within full 3D 
elastodynamics framework [151, 152] suggests promising developments in this context 
in a close future (see [153, 154] for recent theoretical work in this context). 

As presented in section 6.2 this stochastic LEFM description seems to capture 
fairly well the self-afhne scaling features exhibited by fracture surfaces. In particular, 
the apparent disagreements between its predictions and experimental observations 
presented in Bouchaud's 1997 review [10] are now uncovered [18, 95]: The "classical" 
self-afhne scaling regime characterized by a roughness exponent £ ~ 0.8 widely 
reported in the literature is in fact observed below the FPZ size, where in essence 
stochastic LEFM descriptions stop to be relevant. 

While the problem of crack propagation in disordered brittle materials has been 
widely addressed theoretically, the role of disorder on the formation of seed cracks 
remains far less studied (see [155, 156] for recent theoretical work in this topic). 
The fact that fracture surfaces observed at small scale exhibit a universal small scale 
£ ~ 0.8 Family- Viseck scaling in very different materials with various damage processes 
as, e.g., plastic deformation, crack blunting, ductile cavity growth or microcracking 
remains largely unexplained. It suggests that it may be possible to find a generic 
unified statistical description of damage spreading within disordered brittle materials 
within the FPZ, independent of the precise nature of the considered material. In 
this context, damage descriptions in term of universal stress- weighted percolation 
processes as proposed by Schmittbuhl and Hansen (2003) [135] appear promising. 
Their development, their implications in terms of damage dynamics, and their careful 
confrontation to experiments and/or simulations represent interesting challenges for 
future investigations. 
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